Parton energy loss limits and shadowing in Drell-Yan dimuon production 
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A precise measurement of the ratios of the Drell-Yan cross section per nucleon for an 800 GeV/c 
proton beam incident on Be, Fe and W targets is reported. The behavior of the Drell-Yan ratios 
at small target parton momentum fraction is well described by an existing fit to the shadowing 
observed in deep-inelastic scattering. The cross section ratios as a function of the incident parton 
momentum fraction set tight limits on the energy loss of quarks passing through a cold nucleus. 
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The Drell-Yan process, where a beam quark (anti- 
quark) fuses with a target antiquark (quark) producing 
a muon pair, can be used to study the interactions of 
fast partons penetrating through cold nuclei. Only ini- 
tial state interactions are important in Drell-Yan since 
the dimuon in the final state does not interact strongly 
with the partons in the medium. This makes Drell-Yan 
scattering an ideal tool to study energy loss of fast quarks 
in nuclear matter by comparing the observed yields from 
a range of nuclear targets. The dynamics of fast parton 
energy loss in nuclear matter is the subject of consider- 
able theoretical interest jll-^ and has significant implica- 
tions for the physics of relativisitic heavy ion collisions. 

Drell-Yan scattering is closely related to deep-inelastic 
scattering (DIS) of leptons, but unlike DIS it can be 
used specifically to probe antiquark contributions in tar- 
get parton distributions. When DIS on nuclei occurs at 
X < 0.08, where x is the parton momentum fraction, 
the cross section per nucleon decreases with increasing 
nuclear number A due to shadowing pLpl. Shadowing 
should also occur in Drell-Yan dimuon production at 
small X2, the momentum fraction of the target parton, 
and theoretical calculations indicate that shadowing in 
the two reactions has a common origin M,|| . This should 
be particularly apparent at a; < 0.06 where DIS on nuclei, 
like Drell-Yan, is dominated by scattering off sea quarks. 

Fermilab Experiment 866 (E866) measured the nu- 
clear dependence of Drell-Yan dimuon production by 800 
GeV/c protons on Be, Fe and W targets at larger values 



of xi, the momentum fraction of the beam parton, larger 
values of xp (~ xi — X2), and smaller values of X2 than 
reached by the previous experiment, Fermilab E772 [BJ. 
The extended kinematic coverage of E866 significantly 
increases its sensitivity to parton energy loss and shad- 
owing. This Letter reports the results. 

The experiment used the same 3-dipole magnet spec- 
trometer that was described in jlOJ. An 800 GeV/c ex- 
tracted proton beam averaging 3 x 10^^ protons per 20 s 
spill bombarded one of three solid targets or an empty 
target frame. The Be, Fe and W targets were 9.4% - 
19% of an interaction length thick. Their relative thick- 
nesses were chosen carefully to match the background 
rates present in the spectrometer. The targets were lo- 
cated far upstream of the main dipole magnet to opti- 
mize the acceptance at low X2 and large xp- During 
a cycle, 2 beam spills were taken on each target and 1 
spill was collected on an empty location. After passing 
through a target, the remaining beam was intercepted by 
a copper beam dump, which was followed by a thick ab- 
sorber that removed hadrons produced in the target and 
the dump. This ensured that only muons traversed the 
spectrometer's detectors, which consisted of four tracking 
stations and a momentum analyzing magnet. The trigger 
required a pair of triple hodoscope coincidences having 
the topology of a muon pair from the target. Typically 
1400 triggers per second were recorded with an average 
live time of 93%. 

Over 130,000 muon pairs with dimuon mass in the 
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FIG. 1. Ratios of the measured cross section per nucleon for Drell-Yan events versus dimuon mass (in GeV/c ), X2, xf and 
X\. The upper (lower) panels show ratios of Fe/Be (W/Be) from the present experiment as solid circles and Fe/C (W/C) from 
E772 as open circles. The errors are statistical only. The solid curves are the predicted cross section ratios, integrated over the 
hidden variables, from leading-order calculations using EKS98 llll] and MRST WM- 



range 4.0 < M < 8.4 GeV/c^ survived the data cuts. 
Random pairs were subtracted from the data using sim- 
ulated random dimuons constructed by mixing single 
muon tracks that were obtained simultaneously from 
prescaled single muon triggers. The normalization factor 
for the random correction was evaluated by comparing 
to the observed rate of same-charge muon pairs in the 
data. The random rate was 9% of the real events, and it 
occurred primarily at low dimuon mass and high pT- A 
residual 3% background from the empty target was also 
subtracted. The solid targets intercepted the beam at 
nearly the same z location, making differences in dimuon 
acceptance negligible. The overall systematic normaliza- 
tion uncertainty in the cross section ratios reported here 
is 1%. 

The Drell-Yan events obtained by E866 extend over 
the ranges 0.01 < X2 < 0.12 and 0.21 < xi < 0.95, with 
(2:2) — 0.038 and {xi) — 0.46. They also cover the range 
0.13 < Xf < 0.93 and provide good px coverage to 4 
GeV/c. Ratios of the cross section per nucleon for Fe to 
Be and W to Be versus dimuon mass, X2, xp and xi are 
shown in Fig. ^, along with similar results from E772 for 
Fe to C and W to C. The figure shows very good agree- 
ment between the experiments for the cross section ratios 
versus X2- The agreement versus other variables is also 
quite good, given the much smaller {X2) and, hence, in- 
creased shadowing of the present data. Note that the A 
dependence observed here and in E772 implies the change 
in the cross section ratios in Fig. H due to the choice of 
Be versus C is small compared to the effect of the dif- 
ference in {x2)- The reduction in the cross section per 
nucleon on the heavy targets, characteristic of shadow- 
ing, is clearly apparent at small X2- A similar reduction. 



part of which could be related to incident parton energy 
loss, is apparent at large xp and xi. However, it is im- 
portant to recognize that the spectrometer acceptance 
coupled to the intrinsic Drell-Yan cross section leads to 
a strong anti-correlation between X2 and xp for the ob- 
served events. Therefore, the events that show the cross 
section reduction at large xp and xi are in general the 
same events that appear in the shadowing region. 

In order to identify the contributions from shadowing. 
Fig. n^ also shows the predicted cross section ratios, in- 
tegrated over the hidden variables, from leading-order 
Drell-Yan calculations using the code EKS98 [O to- 
gether with the MRST parton distribution functions [Q . 
Essentially identical results are obtained using CTEQ5M 
or CTEQ5L parton distributions ^ instead. EKS98 
describes the ratios fA(x,Q'^)/fp{x,Q'^) of the various 
quark flavors in nucleus A, compared to those in the pro- 
ton. It has been tuned to fit the shadowing observed in 
DIS |14| and E772 while conserving baryon number and 
momentum. 

EKS98 provides an unbiased way to separate the ef- 
fects of shadowing and energy loss in the present data 
because it uses a single shadowing function to describe 
the nuclear dependence of all sea quark distributions at 
its initial scale Qo^ = 2.25 GeV^ and only DIS results are 
used to constrain that function for x ^ 0.1. The shape 
and magnitude of the ratios versus X2 are well reproduced 
by the shadowing predictions. Most of the A dependence 
observed in the ratios versus mass, xp and xi can also be 
explained by shadowing at small X2- As a further test of 
the shadowing parametrization, the events in Fig. m have 
been separated into two sets at the median X2 value. The 
cross section ratios differ by up to 6% when the events 
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FIG. 2. Ratios of the measured Drell-Yan cross section per 
nucleon versus px- Ratios of Fe/Be and W/Be from the 
present experiment are shown as sohd circles, and ratios of 
Fe/C and W/C from E772 are shown as open circles. The 
solid curves are shadowing predictions for the present exper- 
iment from leading-order calculations using EKS98 [|ll| and 
MRST [||. 



at low and high X2 - but at the same mass, xp or xi 
- are observed separately. EKS98 also describes these 
differences well. This is the first experimental demon- 
stration that the shadowing observed in Drell-Yan and 
DIS is quantitatively similar. 

Figure || shows the measured ratios of the Drell-Yan 
cross section per nucleon as a function oi px- Unlike the 
longitudinal variables that have strongly correlated ac- 
ceptances, the Pt acceptance in E866 is nearly indepen- 
dent of the other kinematic variables, so the shadowing 
calculations predict essentially constant cross section ra- 
tios versus pt, as shown by the smooth curves in the fig- 
ure. The data demonstrate a clear pT dependence which 
must have an independent origin. The slight reduction 
in the cross section per nucleon observed for heavy nu- 
clei at small px, coupled with the increase in the cross 
section per nucleon at large pr, is characteristic of mul- 
tiple scattering of the incident partons as they traverse 
the nucleus. 

The xi dependence of the cross section ratios provides 
the best direct measure of the energy loss of the incident 
quarks in the nuclear medium. Table H gives the ratios 
of the measured cross section per nucleon as a function 
of xi, and the mean values of xi, X2 and dimuon mass 
for each bin. However, as shown above, shadowing at 
small X2 explains a substantial fraction of the apparent 
variation in the cross section ratios versus xi. This must 
be removed before one can isolate a nuclear dependence 
due to energy loss. Figure shows the same cross sec- 
tion ratios versus xi as given in Table 0, but corrected 
for shadowing by weighting each event with the calcu- 
lated ratio of the Drell-Yan cross sections per nucleon for 



TABLE I. Ratios of the cross section per nucleon for Fe 
to Be and W to Be as functions of x\ , without correction for 
shadowing. Also given are mean values for xi , X2 and dimuon 
mass (in GeV/c^) for each bin. Mean xf may be found from 
{xf) = {xi) — {x2)- The errors are statistical only. 



x\ range 



(^i) (^2) (Af) 



7-= 



7-' 



0.21-0.3 0.279 0.055 4.76 1.059(41) 0.971(37) 

0.3-0.4 0.356 0.046 4.92 0.973(16) 0.964(15) 

0.4-0.5 0.448 0.039 5.07 0.959(14) 0.951(13) 

0.5-0.6 0.545 0.034 5.17 0.962(17) 0.944(16) 

0.6-0.7 0.642 0.030 5.24 0.979(25) 0.914(23) 

0.7-0.8 0.738 0.026 5.27 0.920(42) 0.884(40) 

0.8-0.95 0.836 0.023 5.28 0.960(99) 0.822(86) 
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FIG. 3. Ratios of the cross section per nucleon versus x\ for 
Fe/Be (upper panel) and W/Be (lower panel), corrected for 
shadowing. The solid curves are the best fit using the energy 
loss form am , and the dashed curves show the Icr upper limits. 
The dotted curves show the la upper limits using the energy 
loss form (H). The la upper limit curves using the energy loss 
form (pl) are essentially identical to those using form (|3|). 



deuterium and nucleus A at the same (xi, X2, Q^), using 
EKS98 and MRST. 

Several groups have studied energy loss of partons in 
nuclei. Their results can be expressed in terms of the 
average change in the incident parton momentum frac- 
tion, Axi, as a function of target nucleus. Gavin and 
Milana M analyzed the E772 Drell-Yan data for energy 
loss based on the parametrization 
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where the factor ki may have a Q^ dependence. They 
based the form (|l|) on an analogy to the transverse spin 
asymmetry in direct photon production. From a com- 
parison to the xp dependence seen by E772 and neglect- 
ing shadowing, they concluded that the fractional energy 
loss of quarks passing through nuclei is « 0.4%/fm. Re- 
cently, other groups have assumed equivalent formula- 
tions with even more rapid energy loss |I3,16|. Brodsky 



and Hoyer [gj argued that the energy loss found by Gavin 
and Milana was too large since the time scale for QCD 
bremsstrahlung was too short to allow for multiple con- 
tributions to the energy loss. Brodsky and Hoyer used an 
analogy to the photon bremsstrahlung process to obtain 
a form for gluon radiation leading to an initial parton 
energy loss [|7| 



Axi 



K2 



A^/\ 



(2) 



They found an upper limit for the gluon radiation from 
the uncertainty relation. They also noted that elastic 
scattering should make a similar contribution to the en- 
ergy loss. Overall, they concluded that energetic partons 
should lose ^ 0.5 GeV/fm in nuclei. The formulation de- 
veloped by Brodsky and Hoyer was extended by Baier et 
al. |^,Q . They found that the energy loss of sufficiently 
energetic partons depends on a characteristic length and 
the broadening of the squared transverse momentum of 
the parton. For finite nuclei, both factors vary as A^/^ ^ 
so Baier et al. predict 



Axi 
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Baier et al. predict that energy loss may be different in 
hot and cold nuclear matter and that there could a large 
coherent effect in relativistic heavy ion collisions. In con- 
trast, tight limits have been placed on energy loss in S-l-S 
and Pb+Pb collisions flj], but the parton momenta rel- 
ative to the hadronic medium were very much smaller 
in those cases than in the present experiment, making a 
direct comparison difficult. 

Given these energy loss expressions, it is possible to 
obtain empirical values for the k's by performing simul- 
taneous fits to the Fe/Be and W/Be Drell-Yan cross 
section ratios versus xi in Fig. ^. We assume that ki 
is Q^-independent ||l|,|[ when using form (|l|) and that 
our incident quarks are sufficiently energetic when us- 
ing form (|^). Curves corresponding to the fits are in- 
cluded in Fig. 0. When assuming the form (nf), we find 
Ki = 0.0004 ±0.0009. This implies that the observed frac- 
tional energy loss of the incident quarks is <0.14%/fm. 
For the energy loss forms (||) and (^, the best fits im- 
ply essentially zero energy loss. We find the Icr upper 
limits to be K2 < 0.75 GeV^ and K3 < 0.10 GeV^. The 
K2 limit indicates that the incident quarks lose energy 
at a constant rate of < 0.44 GeV/fm. The K3 limit im- 
plies that the observed energy loss of the incident quarks 
within the model of Baier et al. is /^E < 0.046 GeV/fm^ 
X L^, where L is the quark propagation length through 
the nucleus. This is very close to the lower value given by 
Baier et al. for cold nuclear matter |l9[ . In all three cases, 
the quoted errors include both statistics and the overall 
normalization uncertainty, with the latter dominating. 

One can also obtain an indirect estimate of the energy 
loss due to gluon radiation in the model of Baier et al. 



from the broadening of the (pt^) as the incident quark 
passes through the nucleus, as shown in Fig. g. How- 
ever, such an analysis involves two significant difficulties. 
While the Drell-Yan cross section is very small for pT 
beyond 4 GeV/c, the change in (pt^) from nucleus to 
nucleus, being the small difference of large numbers, is 
quite sensitive to the yield at very large pr where our 
acceptance is poor and the random background becomes 
significant. Meanwhile, the pt dependence of the ratio 
of the Drell-Yan cross sections per nucleon on hydrogen 
and deuterium shows possible evidence for a change in 
the reaction mechanism at px ~ 3 GeV/c p^, compli- 
cating interpretation of the large pt ratio data in Fig. y. 
Further analysis of the pr dependence of the Drell-Yan 
cross section will be presented in a future publication. 

In summary, this Letter reports a measurement of the 
ratios of the Drell-Yan cross section per nucleon for Fe 
to Be and W to Be. Nuclear shadowing is found to be 
important in the small X2 domain. For sufficiently small 
X2 , the shadowing observed in Drell-Yan has been demon- 
strated to be quantitatively similar to that in DIS. Sub- 
sequently, a correction for shadowing has been applied. 
The cross section ratios versus xi provide direct upper 
limits on the energy loss of the incoming parton as it 
traverses a cold nucleus that are tighter than previous 
constraints. Shadowing and initial state energy loss are 
processes that occur in both Drell-Yan production and 
J/t/j formation. Hence, these results should also further 
the understanding of J/ip production, which is required 
if it is to be used as a signal for the quark-gluon plasma 
in relativistic heavy ion collisions. 
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